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Abstract 

Effective capacity, which provides the maximum constant arrival rate that a given service process can support 
while satisfying statistical delay constraints, is analyzed in a multiuser scenario. In particular, the effective capacity 
region of fading multiple access channels (MAC) in the presence of quality of service (QoS) constraints is studied. 
Perfect channel side information (CSI) is assumed to be available at both the transmitters and the receiver. It is initially 
assumed the transmitters send the information at a fixed power level and hence do not employ power control policies. 
Under this assumption, the performance achieved by superposition coding with successive decoding techniques is 
investigated. It is shown that varying the decoding order with respect to the channel states can significantly increase 
the achievable throughput region. In the two-user case, the optimal decoding strategy is determined for the scenario in 
which the users have the same QoS constraints. The performance of orthogonal transmission strategies is also analyzed. 
It is shown that for certain QoS constraints, time-division multiple-access (TDMA) can achieve better performance 
than superposition coding if fixed successive decoding order is used at the receiver side. 

In the subsequent analysis, power control policies are incorporated into the transmission strategies. The optimal 
power allocation policies for any fixed decoding order over all channel states are identified. For a given variable 
decoding order strategy, the conditions that the optimal power control policies must satisfy are determined, and an 
algorithm that can be used to compute these optimal policies is provided. 



I. Introduction 

In wireless networks, the design and analysis of efficient transmissions strategies have been of significant 
interest for many years. In particular, fading multiple access channels (MAC) have been extensively studied 
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from an information-theoretic point of view [l]-[8]. For instance, Tse and Hanly [4] have characterized the 
capacity region of and determined the optimal resource allocation policies for multiple access fading channels. 
They have shown that the boundary surface points are in general achieved by superposition coding and 
successive decoding techniques, and obtaining each boundary point can be associated with an optimization 
problem in which a weighted sum rate is maximized. Vishawanath et al. [7] derived the explicit optimal 
power and rate allocation schemes (similar to waterfilling) by considering that the users are successively 
decoded in the same order for all channel states. For the convex capacity region, the unique decoding 
order was shown to be the reverse order of the priority weight. While superposition coding and successive 
decoding strategies provide superior performance, time-division multiple access (TDMA) may in certain 
cases be preferred due to its simplicity. Note that the performance of TDMA approaches that of the optimal 
strategy as the signal-to-noise ratio (SNR) vanishes but, as shown by Caire et al. in [8], TDMA is strictly 
suboptimal when SNR is low but nonzero. 

While establishing the fundamental performance limits, the above-mentioned studies have not explicitly 
taken into account buffer constraints and random arrivals. In [9] and [10], Yeh and Cohen considered 
multiaccess fading channels with random packet arrivals to buffered transmitters, and characterized rate and 
power allocation strategies that maximize the stable throughput of the system. In [11], the same authors 
investigated rate allocation policies that minimize the average packet delay in multiaccess fading channels 
again under the assumption of randomly arriving packets. 

In this paper, we also investigate the performance under buffer constraints but provide a perspective 
different from those of previous studies. In particular, we consider statistical quality of service (QoS) 
constraints in the form of limitations on the buffer violation probabilities, and study the achievable rate 
region under such constraints in multiaccess fading channels. Note that in certain delay sensitive applications, 
such as interactive or streaming video, constraints on delay bound violation probability may be required 
rather than limitations on the average delay. For this analysis, we employ the concept of effective capacity 
[12], which can be seen as the maximum constant arrival rate that a given time- varying service process can 
support while satisfying statistical QoS guarantees. Effective capacity formulation uses the large deviations 
theory and incorporates the statistical QoS constraints by capturing the rate of decay of the buffer occupancy 

2 



probability for large queue lengths. The analysis and application of effective capacity in various settings has 
attracted much interest recently (see e.g., [13]— [16] and references therein). In [16], Liu et al. considered a 
two-user cooperative multiple access fading channel and analyzed the rate region achieved with frequency- 
division multiplexing when the users are operating under QoS constraints in the form of limitations on 
buffer overflow probabilities. In this study, cooperation among the users is shown to significantly improve 
the achievable rate region if the quality of the wireless link between the users is better than those of 
the links between the users and the destination. We note that since the transmitters are assumed to not 
know the channel conditions, power and rate adaptation policies are not studied in [16]. Additionally, since 
orthogonal transmission schemes are considered, superposition coding and successive decoding strategies 
are not addressed in detail. 

Our contributions and major findings in this paper can be summarized as follows. We consider the scenario 
in which both the transmitters and the receiver have perfect channel side information (CSI). First, assuming 
that no power control is employed in the transmission, we characterize the rate regions for both superposition 
transmission strategies and TDMA. Unlike the results obtained in [1] and [7], varying the decoding order 
with respect to the channel states is shown to significantly increase the achievable rate region (i.e., throughput 
region) under QoS constraints. Also, it is demonstrated that time sharing strategies among the vertex of the 
rate regions can no longer achieve the boundary surface of the throughput region. Additionally, we show 
that if we take the sum-rate throughput, or the sum effective capacity, as the performance metric, TDMA 
can in certain cases even achieve better performance than superposition coding when fixed decoding order is 
employed at the receiver. Next, we incorporate power control policies into the model. For this case, we first 
obtain closed-form expressions for the optimal power control policies under the assumption that the decoding 
order is fixed at the receiver side. When the decoding order is variable, we identify which conditions the 
optimal power control policies should satisfy. We also describe an algorithm to determine such policies. 

The remainder of the paper is organized as follows. Section II describes the system model. In Section III, 
effective capacity as a measure of the performance under statistical QoS constraints is briefly discussed, and 
the throughput region under QoS constraints is defined. In Section IV, under the assumption of no power 
control, we analyze the throughput region for both fixed and variable decoding order strategies. Section V 

3 




describes the optimal power control policies. Finally, Section VI concludes the paper. 

II. System Model and MAC Capacity Region 

As shown in Figure [H we consider an uplink scenario where M users with individual power and buffer 
constraints (i.e., QoS constraints) communicate with a single receiver. It is assumed that the transmitters 
generate data sequences which are divided into frames of duration T. These data frames are initially stored 
in the buffers before they are transmitted over the wireless channel. The discrete-time signal at the receiver 
in the i th symbol duration is given by 

M 

^^Wl+nl'l 2 = 1,2,... (1) 

j=i 

where M is the number of users, Xj[i) and hj[i) denote the complex- valued channel input and the fading 
coefficient of the jth user, respectively. We assume that {^[«]}'s are jointly stationary and ergodic discrete- 
time processes, and we denote the magnitude- square of the fading coefficients by Zj[i) = \hj[i}\ 2 . Above, n[i) 
is a zero-mean, circularly symmetric, complex Gaussian random variable with variance E{|n[i]| 2 } = N . 
The additive Gaussian noise samples {n[i]} are assumed to form an independent and identically distributed 
(i.i.d.) sequence. Finally, Y[i) denotes the received signal. 

The channel input of user j is subject to an average energy constraint E{|a;j[i]| 2 } < Pj/B for all j, where 
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B is the bandwidth available in the system. This formulation indicates that user j is subject to an average 



power constraint of Pj. With these definitions, the average transmitted signal to noise ratio of user j is 

Mmx/ if \ 

N B- 



p. 

SNRj = jj^. Now, if we denote Pj[i\ as the instantaneous transmit power in the ith frame, the instantaneous 



transmitted SNR level becomes fij[i] = Then, the average power constraint is equivalent to the average 
SNR constraint E{/ij[«]} < sm.j for user j. 

A. Fixed Power and Variable Rate 

First, we consider the case in which the transmitters operate at fixed power and hence do not employ any 
power adaptation policies. The capacity region of this channel is given by [1], [4]: 



72-MAC = | (Ravg,l, Ravg^hl) '■ ^avg(S) < BE Z ^log 2 ^1 + ^ SNR j Z j^j j 



ySc{l,...,M}\ (2) 



where snrj = Pj/(N B) denotes the average transmitted signal-to-noise ratio of user j, z = • • • , z M ) is 
a random vector comprised of the magnitude-squares of the channel coefficients. As well-known, there are 
Ml vertices of the polyhedron defined in ©. The vertex R a ^ j7r = (Ravg,n(i), • • • , Ravg,Tr(M)) corresponds to 
a permutation n, or the successive decoding order at the receiver, i.e., users are decoded in the order given 
by tt(1), • • • , 7r(M). This vertex is specified by the average rates 

Ravg,n(k) = BE Z \ log 2 ( 1 + S ^ (fc) ^ (fc) \ I ^ (3) 

I V 1 + 22i=k+l SNR 7r(i)^7r(i) / J 

With this characterization, we see that for the given decoding order n, the maximum instantaneous service 
rate for user n(k) is 

i^ (fc) = Slog 2 [l + J^W^W ) bits/s k = l,---,M. (4) 

V 1 + Li = fc + l SNR 7T(i)^W/ 

Finally, we note that time sharing among these M\ permutations of decoding orders yields any point 
on the boundary surface of 7^mac [18]. As also discussed in [7], it can be easily verified that varying the 
decoding order according to the channel states does not provide any improvement on the capacity region. 
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B. Variable Power and Variable Rate 

Now, we suppose that dynamic power and rate allocation is performed according to time-variations in 
the channels. For a given set of power allocation policies U = {/ii, • • • , Hm}> where fij > is the power 
control policy of the jth user, the achievable rate region is described by [4] 

K{U) = ^R avg : R avg (S) < E z |filog 2 ^1 + Y^ H (z)z^j J , VS C {1, ■ ■ • ,M}|. (5) 

For a given decoding order at the receiver, the individual average and instantaneous rates of the users can 
be obtained similar to © and ©, respectively, with snr replaced by /i. The capacity region is given by 

^mac = |J K(U) (6) 

where T is the set of all feasible power control policies that satisfy the average power constraint 

7 = {U : E z {^(z)} < SNRj , /ij > 0, Vj} (7) 

where snr^ = Pj/(N B) denotes the average transmitted signal-to-noise ratio of user j. 

C. TDMA 

For simplicity, we assume that the time division strategy is fixed prior to transmission. Let 5j denote the 
fraction of time allocated to user j. Note that we have Y^jLi $j = 1- hi each frame, each user occupies 
the entire bandwidth to transmit the signal in the corresponding fraction of time. Then, the instantaneous 
service rate for user j is given by 

(SNR ■ \ 
1 + -j^-Zjj bits/s - (8) 
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III. Preliminaries 

A. Effective Capacity 

In [12], Wu and Negi defined the effective capacity as the maximum constant arrival ratJJ] that a given 
service process can support in order to guarantee a statistical QoS requirement specified by the QoS exponent 
9. If we define Q as the stationary queue length, then 9 is the decay rate of the tail distribution of the queue 
length Q: 

lim = _ 9 . (9) 

q-+oc q 

Therefore, for large g max , we have the following approximation for the buffer violation probability: P(Q > 
<?max) ~ e _6,l3max . Hence, while larger 9 corresponds to more strict QoS constraints, smaller 9 implies looser 
QoS guarantees. Similarly, if D denotes the steady-state delay experienced in the buffer, then P(D > 
d max ) ~ e-^^max f or laj-ge d max , where £ is determined by the arrival and service processes [14]. Since 
the average arrival rate is equal to the average departure rate when the queue is in steady- state, effective 
capacity can also be seen as the maximum throughput in the presence of such constraints. 
The effective capacity is given by 

C{9) = - lim -^log e E{e- e5[t] } bits/s, (10) 

t-^-oo 9t 

where the expectation is with respect to S[t] = ^2l =1 s[i], which is the time-accumulated service process. 
{s[i],z = 1,2, . . .} denote the discrete-time stationary and ergodic stochastic service process. 

In this paper, in order to simplify the analysis while considering general fading distributions, we assume 
that the fading coefficients stay constant over the frame duration T and vary independently for each frame 
and each user. In this scenario, s[i] = TR[i], where R[i] is the instantaneous service rate in the ith frame 
duration [iT\ (i + 1)T). Then, (TTOb can be written as 

C(9) = -^ f \og e E z {e- 0TR M} bits/s, (11) 

'For time- varying arrival rates, effective capacity specifies the effective bandwidth of the arrival process that can be supported by the channel. 
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where R[i] is in general a function of the fading state z. (fTTT) is obtained using the fact that instantaneous 
rates {-R[i]} vary independently from one frame to another . It is interesting to note that as 9 — > and hence 
QoS constraints relax, effective capacity approaches the ergodic rate, i.e., C(9) — > K z {R[i]}. 

Throughout the rest of the paper, we use the effective capacity normalized by bandwidth B, which is 
denoted by 

C(0) = bits/s/Hz. (12) 

B 

B. Throughput Region 

Suppose that 6 = (9\, ■ ■ ■ ,9 M ) is a vector composed of the QoS constraints of M users. Let C(0) = 
(Ci(6 l i), • ■ ■ , Cm(9m)) denote the vector of the normalized effective capacities. We first have the following 
characterization. 

Definition 1: The effective throughput region is described as 

C MAC (e,SNR) = |J jc(6) > : CjiPt) < -^log e E z {e- 6TR i} J (13) 

s.t. E{R}e^ MAC 

where R = {R\, R2, • • • , Rm} represents the vector composed of the instantaneous transmission (or equiv- 
alently service) rates of M users. Note that the union is over the distributions of the vector R such that the 
expected value E{R} lies in the MAC capacity region. 

Remark 1: The throughput region given in Definition [T] represents the set of all vectors of constant arrival 
rates C(9) that can be supported in the fading multiple access channel in the presence QoS constraints 
specified by 6 = (#1, • • • , 9m)- Since reliable communications is considered, the arrival rates are supported 
by instantaneous service rates whose expected values are in the MAC capacity region. For instance, in the 
absence of power control, the maximum instantaneous service rates for a given decoding order are given by 
©. 

Using the convexity of the MAC capacity region IZmac, we obtain the following preliminary result on 
the effective throughput region defined in (fT3t . 

Theorem 1: The throughput region Cmac(©?snr) is convex. 
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Proof: Let the vectors C(0) and C'(O) belong to C M ac(@,snr). Then, there exist some rate vectors R and 
R' for C(6) and C'(0), respectively, such that E{R} and E{R'} are in the MAC capacity region. By a 
time sharing strategy, for any a G (0, 1), we know from the convexity of the MAC capacity region that 
E{aR + (1 — cc)R'} G IZmac- Now, we can write 

aC(6) + (l-a)C'(e) 

< "Q^g log e (E{e-°™}) a (E{e- e ™'}) 1_a (14) 

1 , ( ' l ^ 1-a 

J- . / _GITVvI3\ — I \ 



QTB 



log e (E{(e- 0TaR ) i }) Q (E|( e - 0T ( 1 - Q ) R ') 1 - Q }) ° (15) 



<-@^B 1 °ge E { e " 0T(ClR+(1 " Q)R ' ) }- 



(16) 



Above, in (fl4l) through (fToT ). all operations, including the logarithm and exponential functions and ex- 
pectations, are component- wise operations. For instance, the expression in (fl4l denotes a vector whose 

components are < log e (E{e"^ T ^}) Q ^E |e _6,Tfi j |^ > . Similarly, the inequalities in (fl4l) and 

L J j=i 

(fToT) are component- wise inequalities. The inequality in (fl4l) follows from the definition in (fT3t . Moreover, 

(fToT) follows from Holder's inequality and leads to the conclusion that aC + (1 — a)C still lies in the 

throughput region, proving the convexity result. □ 

We are interested in the boundary of the region Cmac(@,snr). Now that Cmac(©,snr) is convex, we can 

characterize the boundary surface by considering the following optimization problem [4]: 

max A • C(0) subject to: C(9) G C MA c(0,snr). (17) 

for all priority vectors A = (Ai, • • • , Am) in $H+ J with J2jLi \ = 1- 

IV. Transmissions without Power Control 

In this section, we assume that the signals are transmitted at a constant power level in each frame and 
hence power adaptation with respect to the fading states is not performed. Under this assumption, we initially 
consider the scenario in which the receiver decodes the users in a fixed order. Subsequently, we analyze the 
case of variable decoding order. 



A. Fixed Decoding Order 

We first assume that the receiver decodes the users in a fixed order in each frame. Hence, the decoding 
order does not change with respect to the realizations of the fading coefficients. If a single decoding order is 
used in the frame, it is obvious that only the vertices of the boundary region can be achieved. We consider 
a slightly more general case in which time sharing technique is employed in each frame among different 
decoding orders. Note that the time sharing strategy is also independent of the channel states and hence is 
fixed in different blocks. We denote the fraction of time allocated to decoding order n m as r m . Naturally, the 
fractions of time satisfy r m > and J2m=i T m = 1- Varying the values of r m enables us to characterize the 
throughput region. Under these assumptions, the effective capacity for each user on the boundary surface is 

CM) = -0^ h ^ E * {e-° iT ^ lTmR '~ 1 <fi} (18) 

where R^-i,., represents the maximal instantaneous service rate of user j at a given decoding order 7r m , 
which is given by 

R. 1(i) = B log 2 ( 1 + SNRj ^ ) (19) 

where 7T" 1 is the inverse trace function of 7r m . 

Remark 2: Note that is the maximum instantaneous service rate achieved with superposition 

coding and a particular decoding order. Hence, the corresponding effective capacities characterize the 
throughput achieved with this strategy in the presence of QoS constraints. 

Remark 3: Throughout the rest of the paper, we generally specify the effective capacity values on 
the boundary surface for simplicity and brevity. Effective capacity regions can immediately be specified 
using these boundary points. For instance, the effective capacity (or equivalently throughput) region for 
superposition coding and fixed decoding order is 

|J |c(6) > : CM) < -^log e E z {e-^^-^^O)}] (20) 

{T m } I 3 ) 

where the union is over different time allocation strategies. 
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Next, for comparison, we consider the TDMA case in which we also have similar time allocation strategies 
but only one user transmits in its specific fraction of time. We first have the following definition. 

Definition 2: The throughput region for TDMA can be seen as the achievable vectors of arrival rates with 
each component bounded by the effective capacity obtained when the instantaneous service rate is given by 
©. More specifically, the maximum effective capacity for user j is 

{ { SNR 

1 . „ -<S,-0,-TBlog, ll+^-J-z, 



Cj D (^) = -— log e E<(e ~V" >* "V ) (21, 



'3- 



where 8j is the fraction of time allocated to user j, and < 8j < 1. 
An immediate result can be obtained as follows: 
Theorem 2: The throughput region for TDMA is convex. 

Proof: Note that the points on the boundary surface is given in (|2T1l . Consider the function f(5) = 
—59TB\og 2 (l + ^p- 2 )- It can be easily verified that f(5) is a convex function in 5. Then, e^ s ' is a log- 
convex function. Since weighted non-negative sum preserves the log-convexity [19, Section 3.5], we know 
that E z {e^ 5 )} is log-convex. Then — log e E{e~ <5eT ' Blog2 ( 1H ~ z )} is a concave function in 5. Hence, we 
immediately see that the throughput region for TDMA is convex. □ 

The optimal time allocation policy that maximizes the weighted sum can be obtained through the opti- 
mization problem 

M , ( ( SNR, \ "1 M 
Xj I -SjQjTB log 2 1+ J -- 



^H-Jjk l °^{ e ^ ^ J }- SX - E*J = M,>0. (22) 

The objective function in the above problem is concave, and we can use the Lagrangian maximization 
approach. Taking the derivative of the Lagrangian function with respect to 5j, we obtain the following 
optimality condition for each user: 

SNR, \ / SNR,- 

;nf 



OJ 



E \ e 1 J 7 I log 2 (l + -x^i) snr— !og2 e 

1H~ Jt . ^7 



7vT ~ X i— ? " 7 snr ^ 1 — - K = ( 23 ) 
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where k is the Lagrange multiplier whose value is chosen to satisfy the constraint J2jLi $j — 1- ^ me optimal 
value of 5j turns out to be negative, then the optimal value of 5j should be 0. When Ai = A2 = • • • = \m, 
the obtained values of {Sj} are the ones that achieve the maximal sum-rate throughput, i.e., the sum of 
the effective capacities of the users. Although obtaining closed-form solutions is unlikely, the maximization 
problem in (l22l) can be easily solved numerically using convex optimization tools. Numerical results are 
provided in Section IIV-C1 

B. Variable Decoding Order 

We now study the case in which the receiver varies the decoding order with respect to the fading states 
z = (zi, . . . ,Zm)- More specifically, we assume that the vector space of the possible values for z is 
partitioned into Ml disjoint regions {Z m }^' =1 with respect to decoding orders {n m }^' =1 . Hence, each region 
corresponds to a unique decoding order. For instance, when z £ Z\, the receiver decodes the information 
in the order t\\. Now, for a given partition {Z m }^ =l , the maximum effective capacity that can be achieved 
by the jth user is 

C,(^) = -^glog e E z {e-^} (24) 

= log e ( V / e- e > TR ^ Pz (z)dz) for j = 1,2,..., M (25) 

where p z is the distribution function of z and R^-i^ is given in (THJ1). Akin to the optimization in (fTTt . the 
optimal partition {Z m }^ =1 that maximizes the weighted sum of the effective capacities can be identified by 
solving the following optimization problem: 

M 

max A • C(6) = max \jCj(9j) (26) 

M A • ( m t -e tr \ 

Note that the optimal partition depends on the weight vector A = (Ai, . . . , Am)- By solving a sequence 
of optimization problems for different values of A, we can trace the boundary of the effective throughput 
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region. 

Considering the expression for effective capacity and the optimization problem in d27l) . we note that finding 
closed-form analytical expressions for the optimal partitions of the channel state space seems intractable for 
a general scenario. With this in mind, we consider a simplified case in which all users have the same QoS 
constraint described by 9. This case arises, for instance, if users do not have priorities over others in terms 
of buffer limitations or delay constraints. 

1 ) Two-user MAC: First, we consider the two-user MAC case and suppose that the two users have the 
same QoS exponent 9. Similar to the discussion in [17], finding an optimal decoding order function can be 
reduced to finding a function z 2 = g{zi) in the state space such that users are decoded in the order (1,2) if 
z 2 < g(zi) and users are decoded in the order (2,1) if z 2 > g(zi). Hence, the function g partitions the space 
of the possible values of z — (zi, z 2 ). With this, the optimization problem in (I26T) becomes 

max Aid + (1 - Ai)C 2 (0, g{z x )) (28) 

9 

where Ci{6, g{z\)) and C 2 (9, g{zi)) are expressed as 

1 / POO poo 

C ^9(zx)) = -^logjj o J e- eTB ^ +s ^p z (z h z 2 )dz 2 d Zl 

00 ( A , SNR 121 



-0TB log, 1- 



+ / / e ^ 1+SNR 2 z 2Jp z (zi,z 2 )dz 2 dzi ) (29) 

Jo Jo 

1 / poo pg(zi) 

C 2 (9, g(z 1 )) = log e f ^ I e- eTB ^ 1+s ^p z (z h z 2 )dz 2 dz, 

+ r r e- eTBl < 1+ ^k) Pz { Zl) z 2 )dz 2 dz\. (30) 

Jo Jg(z!) J 

Note that the maximization in (l28l) is over the choice of the function g(zi). Implicitly, g(zi) should always 
be larger than zero as implicitly implied in d29l and d30l) . In cases in which this condition is not satisfied, 
we need to find a function z\ = f(z 2 ) instead, as will be specified below. 

Theorem 3: The optimal decoding order as a function of the fading state z = (zi,z 2 ) for a specific 
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common QoS constraint 6 in the two-user case is characterized by the following functions: 

(1 + SNRi-2i)i^ - 1 



SNR2 

;i + snr 2 z 2 )K~^ - 1 

SNRi 



if if e [1, oo) and 



, if K G [0, 1) 



(31) 
(32) 



where (3 = £p| and K E [0, oo) is a constant that depends on the weight Ai in (l28t and the values of the 
double integrals in (l29t and (l30b . Note that the function used to partition the state space is either g or / 
depending on the value of K. 

Proof: Suppose that the optimal decoding order is specified by the function z 2 = g{z\). We define 



J{g{ Zl )) = AiCiCMOsi)) + (1 " A a )C 2 (M(zi)) 



(33) 



where p(^i) = <?(£i) + s^(^i)- fl^i) is the optimal function, s is any constant, and r](zi) represents arbitrary 
perturbation. A necessary condition that needs to be satisfied is [20] 



d 
ds 



We define the following: 



oo poo 



'i=l I e 

'0 Jg(z 1 ) 



-6TB log 2 (l+SNRi£i) 



p z (z 1 ,z 2 )dz 2 dz 1 



0. 



(34) 



s=0 



00 r9(*i) flTRW (, | SNR r 



./0 



2 2 



2 ^p z (,21, ,22)^2^1 



(35) 



00 /-p(zi) 



JO 



-6TB log 2 (l+SNR 2 z 2 ) 



poo pco 1 SNR 2 z 2 A 

p m (z 1 ,z 2 )dz 2 dz 1 + / / e Slog2 Wi-J^ , z 2 ) dz 2 d Zl 

Jo Jg(zi) 



(36) 
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By noting that ^p! = T](zi), and from d34l> — (l36l> . we can derive 



ds 

-P 



1 + SNR^i 



6TB \ \ 1 + SNR 2 g(zi) 



1-Ax 



'1 + SNR2^(2;i)) p - 



f | SNR 2 C/(2:i) \ ^\ 
\ 1 + SNRiZi / / 



■ Pm{zu g(zi))v(zi)d z i = (37) 



Since the above equation holds for any 77(21), it follows that 



Ai / / SNRi^i 



1 + - (1 + SNRi^i 



\~P 



6TB (pi \\ 1 + sNR 2 g(zi) 

1 - Al j M , ( ^ , SNR 2 s(Zl 



1 + snr^))^ - 1 + =0 (38) 



6TB (f) 2 \ " V 1 + SNR!^! 

which after rearranging and defining K as follows yields 

-P 



n 



+ (l-Ai) 



\ l + 1+SNR 1Z J - U + SNR 2 5-^lJJ 

Obviously, K > 0. Notice that after a simple computation, d39b becomes 

-P 



K. (39) 



1 + SNRi^i 
1 + SNR 2 g{zi) 



K (40) 



which leads to d3TT) after rearranging. Note here that if K < 1, #(21) < for 21 < g s ^ 1 . Then, the 
expressions in (|2"9"1) and (BOl ) are not well-defined. In this case, we denote the optimal function as Z\ = f(z 2 ) 
instead. Following a similar approach as shown in d29l through d40l) yields (1321 . □ 
Remark 4: Above, we have assumed that the users are decoded in the order (1,2) when z 2 < g(zi) (or 
Z\ > f(z 2 ) if K < 1) and decoded in the order (2, 1) when z 2 > g(z{) (or z\ < f(z 2 ) if K < 1). It is 
interesting note that if we switch the decoding orders in the regions (i.e., if users are decoded in the order 
(1,2) when z 2 > g(z\)), exactly the same partition functions as in (I3TT) and (l32l are obtained due to the 
symmetric nature of the problem. Hence, the structure of the optimal functions that partition the space of 
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channel states (21,22) into two non-overlapping regions do not depend on which decoding order is used in 
which region. 

Remark 5: Although the partition does not depend on the choice of the decoding orders in different 
regions, the performance definitely does. Our numerical computations show that the order selected originally 
at the beginning of our discussion (i.e., using the decoding order (1,2) when z 2 < g(zi) or z\ > ffa)) 
provides a larger throughput region than otherwise. This observation leads to an interesting conclusion. Note 
that partition functions g(z\) in (OTT) and /(Z2) in d32l are linear functions of z\ and zi, respectively. When 
K > 1 and 



user 1 is decoded first and user 2 is decoded last. Hence, for instance, when z\ is much larger than z 2 an d 
user 1 is enjoying much better channel conditions, user 1 is decoded first in the presence of interference 
caused by user 2's received signal. User 2, who has less favorable conditions, is decoded subsequently 
without experiencing any interference. Note that such an operation is the opposite of an opportunistic 
behavior and leads to a more fair treatment of users. This is rather insightful since the users are assumed to 
operate under similar QoS limitations (i.e., they have the same QoS exponent 6). Note that if the decoding 
orders are switched, users having favorable channel conditions will be decoded last and hence experience 
no interference. In such a case, there is a bias towards users with better channel conditions, which leads to 
inefficient performance when both users operate under similar buffer constraints. 

Our observations above have led us to propose the following suboptimal decoding order strategy for a 
scenario with more than 2 users. 

2) Suboptimal Decoding Order: In this section, we consider an arbitrary number of users. When all users 
have the same QoS constraint specified by 6, we propose a suboptimal decoding order given by 



Z2 < g(zi) 



(1 + smiZ^K? - 1 



(41) 



SNR 2 



Arr(l) A„-(2) 



... < 




(42) 
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due to the observation that the user with the largest weight A should be decoded last, and the fact that 
the higher the value of z, the less power is needed to achieve a specific effective capacity. Considering 
a two-user example, we, with this choice of the decoding order, can express the points on the boundary 
surface as 

i / p co poo 

Ci(0) = log e J e- 6TB ^ +s ^ Pz ( Zl , z 2 )dz 2 dz x 



/ / Al e - 0TB1 ° S2 ( 1+ ^^)p z (z 1 ^ 2 )dz 2 ^ 1 ) (43) 
o Jo 



Ao2 

r*CO 



0,(0) = -^log e e^^^^z^dz.dz, 

+ r e^ TB1 ° g2 ( 1+ ™^)p z (^, z 2 )dz 2 dz\ (44) 

C. Numerical Results 

We have performed numerical analysis for independent Rayleigh fading channels with E{z} = 1. In Fig. 
El the throughput region of a two-user MAC is plotted for superposition strategies with different decoding 
ordering methods at the receiver, and also for TDMA. In the figure, the solid and dotted curves provide the 
throughput regions achieved by employing optimal and suboptimal variable decoding orders, respectively, 
at the receiver. Note that in the optimal strategy described by the results of Theorem |3j the receiver chooses 
the decoding order according to the channel states such that the weighted sum of effective capacities, i.e., 
summation of log-moment generating functions, is maximized. We see that the suboptimal strategy described 
in Section ITV-B . 21 can achieve almost the same rate region as the optimal strategy, indicating the efficiency of 
this approach. In the same figure, dot-dashed curve provides the throughput region achieved by employing a 
fixed decoding order for all channel states. Here, we observe that the strategy of using a fixed decoding order 
at the receiver is strictly suboptimal even when the users are operating under similar buffer constraints, and 
varying the decoding order with the respect to the channel gains can significantly increase the achievable 
region. Finally, the throughput region of TDMA is given by the dashed curve. We immediately note that 
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Fig. 2. The throughput region of two-user MAC case. SNRi = SNR2 = dB. 8\ — 62 = 0.01. The solid, dotted, dot-dashed, and dashed lines 
represent the regions achieved with optimal variable decoding order, suboptimal variable decoding order, fixed decoding with time sharing, and 
the TDMA respectively. 

TDMA can achieve some points outside of the throughput region attained with fixed decoding order at the 
receiver side. These numerical results show that markedly different strategies may need to be employed 
when systems are operating under buffer constraints. In the absence of such constraints, the performance 
is captured by the ergodic capacity region which cannot be improved by varying the decoding order with 
respect to the channel states [7]. Hence, using a fixed decoding order at the receiver is an optimal strategy 
when there are no QoS constraints. Moreover, TDMA is always suboptimal with respect to the superposition 
schemes regardless of the decoding-order strategy [8]. 

In Fig. |3 sum-rate throughput, i.e. the sum of the effective capacities, is plotted as a function of the 
QoS exponent 9. Here, we note that as 6 increases, the curves of different strategies converge. In particular, 
TDMA performance approaches that of the superposition coding with variable decoding. Hence, orthogonal 
transmission strategies start being efficient in terms of attaining the sum rate under stringent buffer constraints. 
Note that the sum-rate throughput generally decreases with increasing 6, and we conclude from the figure that 
this diminished throughput can be captured by having each user concentrate its power in a certain fraction 
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e 



Fig. 3. The sum-rate throughput as a function of 9. SNRi = 10 dB; SNR 2 = dB. 



of time in the TDMA scheme. We also see that for approximately 6 > 0.006, TDMA starts outperforming 
superposition transmission when a fixed decoding order is employed at the receiver. Such an observation is 
also noted in the discussion of Fig. El In contrast, we observe that as 9 approaches and hence the QoS 
constraints relax, TDMA is the strategy with the worst performance. Note that when the performance metric 
is the ergodic capacity and hence no queueing constraints are considered, this suboptimality of TDMA with 
respect to superposition strategies is well-known (see e.g., [8]). 

We are also interested in the values of parameter K that appear in the functions in Theorem |3]. In Fig. HI 
we plot K as a function of = t^aT- It ls interesting to note that \og e K seems to be linear with respect 

V. Transmissions with Power Control 

In this section, we analyze the case in which the transmitter employs power control policies in the 
transmission. Similarly as before, we initially investigate the scenario in which the decoding order is fixed 
for all channel states. Subsequently, we study variable decoding order schemes. Note that varying the 
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Fig. 4. K vs. ^. SNRi = VQdB. SNR 2 = dB. $1=62 = 0.01. 

decoding order with respect to the channel states, according to the analysis in Section |IVJ has the potential 
to significantly affect the achievable rates. 

A. Power Control Policy for Fixed Decoding Order in All Channel States 

Here, we characterize the optimal power allocation policies when the decoding order is fixed for all 
channel states. Due to the convexity of Cmao there exist Lagrange multipliers k = (/%,..., km) £ such 
that C*(6) on the boundary surface can be obtained by solving the optimization problem 

max A- C(9) -K-E{fi} (45) 

where fi = . . . , fj, M ) represents the collection of the power control policies of all users, A = (Ai, . . . , Am) 
is the weight vector, and C(0) = (Ci(0i), . . . , Cm (9m)) 1S tne vector of maximum effective capacities of 
the users for given decoding order and power allocation policies. Note that fXj = (defined in Section ITTI 
as the instantaneous transmitted SNR level) describes the power control policy of the jth user . For a given 
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permutation n and set of power allocations \i, Cj(9j) is given by 



,-TB log, 1 



= -_^_log e E^ e 3 V ^-i^-i^wj v, . (46) 



Now, the optimization problem (143b can rewritten as 

max^-A,— -log e E^e J V ^-^-i u) ^J I _ J^E^}. (47) 

The following result identifies the optimal power adaptation policies that solve the above optimization 
problem. 

Theorem 4: Assume that the receiver, for all channel states, decodes the users in a fixed order specified 
by the permutation n. Then, the optimal power allocation allocation policies that solve the optimization 
problem in (|4*71 are given by 

Vi=\- fr L () 0) forj = 1,2,. ..,M (48) 



/ (} + E^-i(p>^-i(j) lazi)"' ^ 1 + E^-i(i)>^-i(j)/^^A 

I _l b_ z I 



TB 

where (3j = is the normalized QoS exponent, = max{i, 0}, and (a 1; • • • , a M ) are constants that 
are introduced to satisfy the average power constraints. 

Proof: Note that with a fixed decoding order, the user tt(M) sees no interference from the other users, and 
hence the derivative of d47l) with respect to ^(m) will only be related to the effective capacity formulation 
of user 7r(M). Therefore, we can solve an equivalent problem by maximizing C^iM) instead. After we derive 
/i^M), the derivative of d47t with respect to /^(Af-i) will only be related to the effective capacity formulation 
of user 7r(M — 1). By repeated application of this procedure, for given A, (|4*71 can be further decomposed 
into the following M sequential optimization problems 



1 -0,TBlog, ^ 

max -A,— — log e E < e \ 1+ ^- 1 (o>-- 1 o) } - /,,?:{//,} j e {1, - • • , M} (49) 
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in the inverse order of ir. Similarly as in [13], due to the monotonicity of the logarithm, solving the above 
M optimizations is the same as solving 



e V 1+ Vi«>»-i(j)«'*/ K + KjE{nj} je{l,---,M}. (50) 

Differentiating the above Lagrangian with respect to fij and setting the derivative to zero yield the intended 
result in d48t . □ 

Remark 6: Exploiting the result in (BHl . we can find that instead of adapting the power according to only 
its channel state as in [13] where a single-user scenario is studied, the user adapts the power with respect 
to its channel state normalized by the observed interference and the noise. 

Remark 7: To give an explicit idea of the power control policy, we consider a two-user example in which 
the decoding order is (2, 1). For this case, we can easily find that 

77 z l > Oil, 

otherwise 

and 



A*i=< ' , (51) 



1 z\ < a\ and z 2 > a 2 , 



- 1 J2 2 2 

in + l lio + 1 

I (±L) (/3i+1)(/3 2 + 1) (ll)WFr , \t±tt 

otherwise 
where a x and a 2 are chosen to satisfy the average power constraints of the two users. 

B. Power Control Policy for Variable Decoding Order 

In this section, we study the optimal power allocation policy when the receiver varies the decoding order 
with respect to the channel fading states. We mainly concentrate on the two-user scenario. The key idea 
we introduce here is to consider the power allocation policy of each user j for each region Z m (in which 
decoding is performed according to permutation n m ) while requiring the average power constraint to be 
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satisfied by the joint power over all regions {Z m }. 

For the two-user case, due to the convexity of the throughput region, there exist Lagrange multipliers 
k = (k,i,k 2 ) G d\\ such that C*(0) on the boundary surface can be obtained by solving the optimization 
problem 

max AiCi(/x, Z) + A 2 C 2 (/i, Z) - KiE{/ii} - K 2 E{/i 2 } (53) 

where /i = (/i 1; /i 2 ) are the power control policies, (Ai,A 2 ) are the weights in the weighted sum, and 
Z = (Zi,Z 2 ) denotes a particular partition of the space of the positive values of z = (zi, z 2 ) Q Hence, 
power control policies that solve (1531 are the optimal ones for a given partition. In the following, since we 
assume Z is given, the notation Cj(fi,Z) is replaced by C,(/z) for brevity. 

Recalling the discussion in Section IIV-B1 we can express the effective capacities of the two users as in 
d29t and (l30l) by only replacing snr, with /i,(z) in these expressions. The Lagrangian (which is the objective 
function in (l53l ) can now be expressed as 

J = ~ q i 1 log e ( / ( 1 + ^ lZl ) p z (z 1 ,z 2 )dz 1 dz 2 + [ (1 + HiZx)' 131 p z (zi, z 2 )dz x dz 2 ) 
Pilog e 2 \J xeZl \ I + J Z £Z 2 J 

\ ( r / \ -&> 

A 2 1 1 / I 1 , A 4222 



log e / I 1 + - — ■ Pz(zi, z 2 )dzidz 2 + / (1 + ^2^2) V*{z\, z 2 )dz l dz 2 



Jz£Z 2 \ 1 + ^1^1/ JzG^i 

- «i(E«e2i{^i} + E ze22 {yUi}) - n 2 (E zeZl {fi 2 } + E zeZ , 2 {fi 2 }). (54) 

Above, the expressions in regions Z\ and Z 2 are written separately due to the reason that possibly different 
power allocation strategies are employed in different regions. We define 

01=/ ( 1 + z ^ lZl ) p z {z 1 ,z 2 )dz 1 dz 2 + (1 + /ii^i)^ 1 p z (zx, z 2 )dzidz 2 , (55) 
J.eZi \ 1 + ^2^2/ J ze z 2 

2 Similarly as discussed in Section HV-BI different decoding orders are employed in Z\ and Z 2 . 
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and 

2 = / ( 1 + ^ 2 * 2 J p*(zi, z 2 )dz 1 dz 2 + / (1 + /i 2 ^)~ /32 p«(zi, z^dsitfefc- (56) 
Jzez 2 \ 1 + fMiZiJ J zeZl 

Note that the values of these functions are obtained for given power control policies /i — (jii, /i 2 ) and given 
partition Z = (Zi, Z 2 ). 

Now, we consider the power control policy of each user in each decoding order region Z^ i = 1,2. By 
differentiating the Lagrangian, we can find the following optimality conditions: 

1 — (1 +/*i*i) *i _ Xl o 1 + 7T 77T ^ _K i = Vze2 i ( 57 ) 

2 log e 2 V I + (l + zii^) 2 

2) T ^(1 + T ^-) -^-« 2 = VzeZ, (58) 





01 log 


; e 2 


A 2 




02 log 


; e 2 


Ai 




01 log 


; e 2 




Ai 


01 


log, 



3) ( 1 1 ^ ) * 1 Kl = Vz G Z 2 (59) 

1 + /i 2 ^ 2 

\ -A-i \ 

^1^1 \ , ^2 n \-/3 2 -l 



+ , , n (l + to^Y* 2 - k 2 = VzG2 2 (60) 



>i log e 2 \ 1 + /i 2 ^ 2 / (1 + /i 2 z 2 ) 2 2 log e 2 

where (1571 and (I5S1) are obtained by differentiating the Lagrangian with respect to \i\ and /i 2 , respectively, 
over z G Z\. Similarly, (1591 and (l6T)t are obtained by differentiating with respect to [L\ and /i 2 , respectively, 
over z G iT 2 . Due to the convexity, whenever n^i = 1,2 is negative valued, we set m = 0, i = 1,2. 
Although obtaining closed form expressions from the optimality conditions seems to be unlikely, we can 
gather several insights on the power control policies by analyzing the equations (l57l-(l6T)t. 

Let us first define e*i = Kl ^ 2 , a 2 = K2 ^ 2 , « 12 = K2 ^ 2 , and a 21 = where Kl , k 2 are 

the Lagrange multipliers whose values are chosen to satisfy the average power constraint © with equality, 
and 0i and 2 are defined in (l55t and (l56b . Now, consider (l57t and (l58t . The channel state lies in Z\. 
Through a simple computation using (l58l) . we can derive 

(l + fi lZl )^ 1 + fj, lZl 
[12 = ^_ ; (61) 

Vt 2 £ 2 
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which tells us that /x 2 = if 



' 2 < a 2 . (62) 



1 + /IiZi 



If fj, 2 = 0, we have from d57t that 



ii log e 2 



1 + fi 1 z i y^ L z 1 - ki = (63) 



which gives us that 



which implies that /ii = if 



Mi = — - y (64) 



0i < ai. (65) 



Now, if we substitute (I6TT) into d57t . we obtain the following additional condition for having [L\ = 0: the 
equation 



£i (1 + _ if*. * - 1 | - l = u (66 ) 

on z 2 a 12 ^"2(1 +ViZi)J 

has a solution that returns a negative or zero value for /ii. The above discussion enables us to characterize 
the regions in which one user transmits while the other one is silent. We also have a closed-form formula 
in (l64l for the optimal power adaptation policy when only one user transmits. Indeed, this is the optimal 
power control policy derived in [13] for a single-user system. When both users transmit, the power control 
policies (/ii,/i 2 ) are given directly by the non-negative solution of (l57l) and (l58t . 

Note that the conditions and characterizations provided in (l6TT) - (l66l) pertain to the case in which the 
channel state is in region Z\. Following a similar analysis of (l59t and (l60l) . we can obtain similar results 
for the cases in which the channel state is in Z 2 . 

For a given partition {Z±, Z 2 }, the optimal power control policy can be determined numerically using the 
optimality conditions in (BTt - (l60l) . Additionally, the equations and inequalities in (I6TT) through d66t can be 
used to guide the numerical algorithms as they specify under which conditions at most one user transmits, 
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and provide the optimal power control policy in such cases. However, there is one difficulty. (1611 - 4661 
depend on a%, a 2 , «i2, and a%\ which in turn depend on <f>i, </> 2 , «i, and k 2 which are in general functions of 
the power control policies. In such a situation, the following iterative procedure can be employed in search 
of the solution . We can first choose certain values for 0i, </> 2 , Ki, and k 2 , and then determine the optimal 
power allocation policies for these selected values. Subsequently, we can check whether the obtained policy 
satisfies the average power constraint with equality. This enables us to determine if the selected Ki and k 2 
values are accurate. We can also compute <pi and 2 using the obtained policy and see if they agree with 
the initial values of <pi and 2 . If there is no sufficient match or if the power constraint is not satisfied with 
equality, then we update the values of 0i, 2 , K\, and k 2 , and reiterate the search of the optimal policy. 

With this insight, we propose the following algorithm that can be used to determine the optimal power 
allocated to each channel state: 
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Power Control Algorithm 

1 Given Ai, A2, the partition Z, initialize 0i, 4> 2 ; 

2 Initialize K\ and k 2 ; 



3 Determine a x = K ^^ 2 , 

4 if z 6 Z x 

5 then if z 2 > a 2 



then /i 2 



Q: 2 



1 



log e 2 
A 2 



"12 



K2</>1 log e 2 



"21 



Kitfe log e 2 . 
A 2 



1 /3?~ 
*2 2 



J_. 

22' 



7 
8 
9 

10 
11 
12 



if f (1 + mi*) 
then /ix = 0; 



else if ^ < Ml 

a 2 \ ai 



-(/3i+l) _ z ia2 
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22«21 \ V a2(l+MlZl) 



02 + 1 



1—1 = returns nonpositive /Xi 



?i+i 



1 



— r~ — /?r~ 



then fi 2 = 0, /ii 
else Compute //i,A*2 from (BTt and d58b : 
else ^2 = 0, /ii = 



1 



— — w 

*1 2 1 



J_ 

21 



13 if z e Z 2 

14 then if z x > «i 



15 

16 
17 
18 

19 
20 
21 



then A^i 



1 

~T- — W~ 



21' 



if 5 (1 + A^ 2 ) 
then /i2 = 0; 



-(/32 + 1) 22Q1 
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zicm \ V ai(l+/i2Z2) 



01+1 



returns nonpositive /i 2 



else if — < ( — 



02 + 1 



then Aii = 0, /i 2 



1 



1 j8 2 
*2 z 2 



J_ 

Z2 



else Compute Aii>A*2 from d59l and d6H: 



else a*i = 0, \i 2 



1 



1 J2 
l 2 z 2 



Z2 



22 Check if the obtained power control policies a«i and fj, 2 satisfy the power constraint with equality; 

23 if not satisfied with equality 

24 then update the values of Ki and k 2 and return to Step 

25 else move to Step EE 

26 Evaluate 0i and 2 with the obtained power control policies; 

27 Check if the new values of </>i and </> 2 agree (up to a certain margin) with those used in Step 

28 if do not agree 

29 then update the values of <\>\ and <p 2 and return to Step El 

30 else declare the obtained power allocation policies [i\ and /i 2 as the optimal ones. 
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Note that we above have not specified how the values of «i, k 2 , 0i, and 02 are updated for each iteration 
in order to keep the algorithm generic. In our numerical computations, we have updated K\ and k 2 using the 
bisection search algorithm. The values of 0i and 02 are updated in Step 29 of the algorithm by assigning 
them the values evaluated in Step 26. Hence, the most recent values are carried over to the new iteration. 

In Fig. we plot the optimal power allocation policies [i\ and /i 2 as functions of channel fading states 
Z\ and z 2 . We assume that 61 = 9 2 = 0.01, snr x = snr 2 = dB, and Ai = A 2 = 0.5. We consider the 
partition specified by the suboptimal decoding order given in (l4"2l . Hence, since we have Ai = A2 = 0.5, 
decoding orders (1,2) and (2,1) are used when z 2 < Z\ and z 2 > Z\, respectively. Under these assumptions, 
we computed the optimal values as k\ = 0.0470, k* 2 = 0.0462, <p\ = 0.5550, and (f> 2 = 0.5538. In the figure, 
we observe that each user, not surprisingly, allocates most of its power to the regions in which it is decoded 
last and hence does not experience interference. However, due to the introduction of QoS constraints, we also 
note that each user also allocates certain power to the cases in which it is decoded first. This is performed 
in order to continue transmission and avoid buffer overflows. 

So far, we have assumed that the partition Z is given. The optimal partition Z that maximizes the weighted 
sum-rate can be derived through the following optimization similarly as in [21]: 

C* = supA 1 C 1 (/i,Z) + A 2 C 2 (/i,Z) (67) 

z 

where C* is the optimal weighted sum value for given pair of (Ai,A 2 ), and /i = (/xi,/i 2 ) are the optimal 
power control policies for given Z. 

VI. Conclusion 

In this paper, we have studied the achievable throughput regions in multiple access fading channels 
when users operate under QoS constraints. We have assumed that both the transmitters and the receiver 
have perfect CSI. We have employed the effective capacity as a measure of the throughput under buffer 
constraints. We have defined the effective capacity region and shown its convexity. We have considered 
different transmission and reception scenarios e.g, superposition coding, different strategies for the decoding 
order, and TDMA. Under the assumption that no power control is employed by the transmitters, we have 
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Fig. 5. The optimal power control policies /ii and fii of users 1 and 2, respectively, as a function of (2:1,2:2). Ai = 0.5, A2 = 0.5. 



29 



analyzed the performances of fixed and variable decoding order strategies. We have characterized the 
throughput region and determined the points on its boundary for fixed decoding order. For the case of 
two users with the same QoS constraints, we have derived the optimal strategy for varying the decoding 
order. Varying the decoding order is shown to significantly increase the achievable rate region. We have 
also proposed a simpler suboptimal decoding rule which can almost perfectly match the optimal throughput 
region. We have also studied the performance of orthogonal transmission strategies by considering TDMA. In 
the numerical results, we have demonstrated that TDMA can perform better than superposition coding with 
fixed decoding order for certain QoS constraints. More specifically, we have noted that TDMA can support 
arrival rate pairs that are strictly outside the region achieved when fixed decoding order is employed at the 
receiver. We have also observed that the performance of TDMA approaches that of the optimal strategy of 
superposition coding with variable decoding order as 9 increases (i.e., as the QoS constraints become more 
stringent). In the second part of the paper, we have incorporated power adaptation strategies into the model. 
For a given fixed decoding order at the receiver, we have identified the optimal power control policies. 
For cases in which a variable decoding order strategy is adopted by the transmitter, we have obtained the 
conditions that the optimal strategies should satisfy and described an algorithm to achieve these optimal 
schemes. 
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